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A simple argument shows that eigenstates of a classically 
ergodic system are individually ergodic on coarse-grained 
scales. This has implications for the quantization ambiguity 
in ergodic systems: the difference between alternative quan- 
tizations is suppressed compared with the Oih?) ambiguity 
in the integrable case. For two-dimensional ergodic systems 
in the high-energy regime, individual eigenstates are indepen- 
dent of the choice of quantization procedure, in contrast with 
the regular case, where even the ordering of eigenlevels is am- 
biguous. Surprisingly, semiclassical methods are shown to be 
much more precise for chaotic than for integrable systems. 

For many years, it has been widely recognized that 
"quantizing" a given classical system is inherently an am- 
biguous procedure, as a large family of quantum Hamil- 
tonians may have the same classical limit [P. For ex- 
ample, given the classical dynamics of a particle con- 
strained to move on a closed loop, different choices of 
boundary condition give rise to different phases relating 
classical paths of different winding number. Knowledge 
of these Aharonov-Bohm phases is of course necessary to 
construct a semiclassical dynamics (which includes inter- 
ference between classical paths) , and thus many semiclas- 
sical theories correspond to the same classical dynamics. 
Physically, this 0{h) or gauge ambiguity may be asso- 
ciated with the possibility of varying the magnetic flux 
enclosed by the loop. 

This is not all, however: there are also many quantum 
theories, differing at 0{fi ) or higher in the Hamiltonian, 
which all have the same semi-classical limit. There are 
many ways of seeing this 0{h ) ambiguity; one of the 
simplest is to imagine making a canonical transforma- 
tion on the classical phase space, applying the canoni- 
cal quantization prescription to the new coordinates, and 
then transforming back to the original coordinate system. 
Generically, one then obtains a new quantum Hamilto- 
nian which differs from the original by 0(?i ) plus higher 
order terms: 



H' =H + h^A 



(1) 



where the operator A has a well-defined classical limit [g . 
As classical dynamics is of course independent of the 
choice of coordinate system, this implies that quantiza- 
tion is inherently ambiguous at second order in h. 

An even more striking case is that of a particle con- 
strained to move on a two-dimensional surface embedded 
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in three-dimensional space. Here, the non-trivial met- 
ric contained in the kinetic term gives rise to obvious 
operator-ordering ambiguities in canonical quantization; 
this has led to much discussion in the literature over 
whether a term proportional to the local Gaussian cur- 
vature R of the surface should be added to the Hamilto- 
nian, and if so, what the proportionality constant should 
be (different prescriptions suggesting Ti R/8, fi R/6, and 
h^R/ 12 as the "correct" answer) ^^. In the path in- 
tegral approach, ambiguities at the same order arise in 
choosing how to incorporate the metric into the kernel 
and in deciding at what point in the infinitesimal time in- 
terval to evaluate functions of the metric . Physically, 
one may define dynamics on a constraint surface through 
a limiting process, where the strength of restoring forces 
causing the particle to live on the surface is taken to in- 
finity M . Classically, this procedure is known to give an 
unambiguous constrained dynamics ||^; in quantum me- 
chanics, 0{h ) differences arise depending on the precise 
way in which the strength of the constraining potential 
is taken to infinity at various places along the surface [pi . 
The ambiguity here has a clear physical meaning: to de- 
termine the true quantum mechanics one needs to know 
the mechanism through which the particle is bound to 
the surface of constraint; it is not sufficient to know only 
the intrinsic properties of the surface itself. Similarly, 
there is no "correct" answer to the problem of quantiz- 
ing a classical double pendulum H]: quantum dynamics 
at 0{h ) is determined by the precise way in which one 
takes to infinity the rigidity of the two rods. 

In a two-dimensional system, the energy spacing be- 
tween adjacent levels is 0{h^), i.e. of the same order 
as the quantization ambiguity demonstrated above. It 
would then seem not to be possible to uniquely determine 
the eigenvalues and eigenstates of a two-dimensional 
system given only the classical dynamics on the two- 
dimensional surface. Similarly, it should not be generally 
possible to compute semiclassically the levels and wave- 
functions of a two-dimensional quantum system, because 
a given semiclassical calculation has many quantum the- 
ories corresponding to it, with Hamiltonians related as 
in Eq. |l|. Recently, however, it was shown (and numer- 
ically confirmed) that in strongly chaotic systems long- 
time semiclassical methods can in fact be used to com- 
pute quantum properties to accuracy much better than 
a level spacing, i.e. that individual wavefunctions and 
eigenenergies can be semiclassically resolved [13]. This 
creates an apparent paradox, to be resolved in the present 
paper. 

The answer is somewhat surprising: it is found that (i) 



the quantization ambiguity is greatly reduced for clas- 
sically ergodic as compared with regular systems; (ii) 
in two dimensions, the ambiguity for ergodic systems is 
small compared to a level spacing, in contrast with the 
integrable case where even the ordering of eigenlevels is 
ambiguous; and (iii) unexpectedly, semiclassical meth- 
ods are valid to much longer times in strongly chaotic 
as compared with integrable systems, allowing individ- 
ual eigenenergies to be easily resolved. The key to these 
surprising results is the coarse-grained ergodicity of indi- 
vidual wavefunctions in classically ergodic systems. More 
precisely, in the classical limit, the quantum expectation 
value of any classically defined operator over individual 
eigenstates must approach the ergodic, microcanonical 
average of the operator, for almost all eigenstates. In 
effect, the fraction of eigenstates that deviate from er- 
godicity when smoothed over a finite mesh size in phase 
space must tend to zero in the ?i — > limit. This be- 
havior has been studied mathematically under a variety 
of technical assumptions by Shnirelman, Zelditch, and 
Colin de Verdiere [Q . Here we present a simple physical 
argument that demonstrates the generality of the result. 
Let A be an arbitrary quantum operator having the 
smooth phase space function A{q,p) as its classical limit. 
Now we select any energy Eq and choose a classically 
small energy window AE such that the microcanonical 
average a{E) of the function A{q,p) is as close as one 
likes to a constant, qq, for all E e [Eo,Eq + AE]. No- 
tice that in the h —^ limit, the classically infinitesimal 
energy window will contain an arbitrarily large number 
of quantum eigenstates, N ^ h~ in d dimensions. We 
wish to show that for any e, there exists h such that 



({n\A\n) 



ao 



<e" 



(2) 



where the average is taken over all wavefunctions \n) in 
the energy window [Eq, Eq + AE]. 
Consider the quantum correlator 



fit) = tr iU(t) = tr A^'e'^'Ae 



Ht A„-iHt 



(3) 



where the trace is over all states in the energy window. 
Now define the time-averaged correlator: 



FiT) 



1 



V2nT 

N 



dte-'^'^^'m 



^ E I (^^1^1"') 



-(£;„-B„/)^T^/2 



(4) 



Fci(T), the time-averaged correlator of A{q,p) and 
A{q{t) , p{t)) , is the classical counterpart of F{T), and 
must by the definition of ergodicity tend to the ergodic 
value Oq as r — > oo. Now we simply choose T large 
enough so that Fc\{T) is as close as we like to its long-time 
asymptotic value (say, within 0(e^)), and then choose 
fi small enough so that the Ehrenfest time Teiii- (the 



time at which classical-quantum correspondence breaks 
down) is large compared with T. [Notice that for a hard 
chaotic system, the breakdown time lEhr ^ A^^log?i~ 
as 7i — > 0, where A is the Lyapunov exponent; however 
we make no specific assumption here about the system 
apart from its being ergodic] We then have 



1 



N 



-^\{n\A\n) <F{^)<F{T) = ai + 0{e'). (5) 

n=l 

Now cq is of course the mean value of the matrix ele- 
ments (n|A|n) (which can be seen formally by consider- 
ing the correlation of A{t) with the identity operator), 
so Eq. ra implies that the variance of the matrix elements 
is bounded by e^ (proving the statement of Eq. ^, and 
goes to zero in the h ^> Q limit. 

In the presence of symmetry-induced degeneracies, the 
argument above clearly holds for any choice of orthonor- 
mal eigenbasis \n). Furthermore, the difference between 
the first two quantities in Eq. || is also bounded above by 
O(e^); this implies 



{n\A\n') 







(6) 
and In') in the 



for almost all degenerate eigenstates 
?7, — > limit. 

As an aside, we mention that the constraint F{oo) < 
F{T) (obtained from Eq. Q) allows us to circumvent the 
usual problem of noncommutativity of the ?i — > and 
T —> oo limits. Specifically, quantum approach to ergod- 
icity at short times in the sense of F{T) -^ Gq guaran- 
tees quantum ergodicity at longer times where individual 
eigenstates are resolved. 

The statistical argument presented here leaves open 
the possibility of a zero measure set of states grossly vi- 
olating the ergodic condition (n|A|n) « gq. This is not a 
deficiency in the argument: infinite sequences of highly 
non-ergodic states (namely, the bouncing-ball wavefunc- 
tions) have been shown to exist in chaotic systems with 
marginally stable classical orbits |I^. The fraction of 
such states does go to zero in the classical high-energy 
limit, in accordance with the predictions of the theorem. 

We also note that no assumptions have been made in 
the argument about "hard chaos" properties such as hy- 
perbolicity or mixing; implications for the elimination of 
quantization ambiguity in two dimensions thus hold for 
all classically ergodic systems. 

Coarse-grained quantum ergodicity is a much weaker 
and more general result than random matrix theory 
(RMT), which has been conjectured to be a statistical 
description of classically ergodic wavefunctions ||l^ . The 
latter conjecture suggests ergodic behavior of quantum 
wavefunctions on single-wavelength scales, and is vio- 
lated by scars |lj] and other effects of short-time dynam- 
ics [nSllGl. For example, the quantum wavefunctions of 
the Sinai billiard, a paradigm of classical chaos, have been 
shown to have inverse participation ratios (IPR's) that 
diverge in the classical limit away from their ergodic and 



RA/[T values |1^. This strongly non-ergodic wavefunc- 
tion behavior on scales of size h is entirely consitent with 
ergodicity on coarse-grained scales as discussed above. It 
is easy to see that the ergodicity argument breaks down 
if A does not have a well-defined classical limit. We may 
for example let v4 be a projection onto a localized test 
state A = \a){a\, but then we are not guaranteed to find 
a time T where classical-quantum correspondence be- 
tween F{T) and Fc\{T) still holds, and simultaneously 

T « TEhr. 

Comparing the result of Eq. with Eq. |l|, we see that 
the leading effect of a change in the quantization pre- 
scription is to shift all energy levels classically close to 
Eq by a constant displacement fi^ao (comparable in size 
to a level spacing); perturbation theory then gives: 



SE,, = E:^-En = h^iao + Oie)) 



(7) 



with e ^ 1 for small h (high energy). In the special 
case of RMT, the energy hypersurface has M ^ h ~ 
Planck-sized cells: since the deviations e in the matrix 
elements (n|v4|n) arise from uncorrelated fluctuations in 
these cells, we then have 



1 



ermt 



-(d-l)/2 



M 



The relative shift between nearby energ; 
to a level spacing is 



(8) 
;y levels compared 



6E„ - 5E„ 



RMT 



^(3-d)/2 



(9) 



The overall energy shift h qq is of course unphysical, 
corresponding simply to changing the Hamiltonian by a 
constant; only energy differences are measurable, and 
from the first relation in Eq. ^ we see that for d = 2 
these become independent of one's choice of quantization 
in the ?i ^ limit (i.e. for highly excited eigenstates). 
Energy splittings £"„ — £",„ which are already classically 
large can of course change by 0{h ) as one considers 
different quantum systems with the same semiclassical 
limit. However, the ratio S{En — Em)/{En — E„i) -^ 
for almost all states \n), \in) in the semiclassical limit for 
d — 2, whether \n) and \m) are nearest neighbor levels or 
are far apart in energy. Thus, in the classical limit of a 
d = 2 constrained surface, the way in which the particle 
is bound to the surface has no effect on any measurable 
quantity, provided that motion on the constraint surface 
itself is ergodic. [In the case of mixed phase space, states 
living on regular islands move up or down relative to the 
rigid chaotic sea, as the quantization is varied.] 

The result of Eq. |9| is consistent with the finding in 
pO| that in caustic-free chaotic d = 2 systems, long- 
time semiclassical dynamics approaches the quantum an- 
swer, and that the critical dimension for breakdown of 
the semiclassical approximation at the Heisenberg time 
is d = 3 for chaotic systems, as compared with d = 2 in 
the integrable case. 



For a simple numerical example of the above results, 
we consider a discrete-time map on a two-dimensional 
toroidal phase space {q,p) G [0,1) x [0,1) (notice that 
this scales equivalently to a two-dimensional autonomous 
Hamiltonian system and can be thought of as a Poincare 
section of the latter) : 



p ^ p = p — V (q) mod 1 
q —>■ q = q + p mod 1 . 



(10) 



The above equations of motion can be obtained from the 
time-periodic ("kicked") Hamiltonian 



H 



2 



V{q)Y.6{t- 



We let the potential be given by 
"2 0.4 



V{q)^±- 



(27r) 



sin 2'Kq + rh cos 2iTq , 



(11) 



(12) 



where the — sign gives completely chaotic dynamics [|0| , 
while the -I- sign leads to a (mostly) regular classical 
phase space. The sin27rg term is added to break parity 
symmetries and make the dynamics nonlinear and generic 
in both cases, while r parametrizes a one-parameter fam- 
ily of possible quantizations. 

A sample piece of the regular (solid curve) and chaotic 
(dashed curve) spectrum is shown in Fig. 1, as a function 
of the quantization parameter r. We see in the regular 
case that the eigenvalues often cross each other as the 
quantization is varied, making impossible a semiclassical 
ordering of the spectrum, while in the chaotic case the 
eigenvalues never shift by an amount comparable to a 
mean level spacing [k, 0.0117 in this calculation). We 
notice that in our example the trace oi A{q,p) = r cos 2'Kq 
is zero, so the overall spectral shift h ao of Eq. vanishes 
as well. For a different A{q,p) = r [cos 27rg ± c] , there 
would be a secular upward or downward trend in all the 
eigenvalues (regular and chaotic) as r increased, but no 
physical quantities would be affected. 

Despite the exponential proliferation of classical paths 
in the chaotic case, semiclassical calculations past the 
Heisenberg time can be performed to arbitrary precision 
with only a power-law amount of effort, using an iter- 
ative approach, and individual semiclassical eigenvalues 
and eigenstates can be extracted [|lO| . Semiclassically ob- 
tained eigenvalues for the chaotic system are indicated by 
arrows in Fig. 1 and are observed to agree well with the 
quantum results (independent of quantization). In the 
regular case, there is of course little meaning to semiclas- 
sically computed individual eigenvalues, and these are 
not shown. 

In Fig. 2 we show the mean squared change in eigen- 
level position, as a function of effective Planck's con- 
stant (inverse momentum) h = 1/64 • • • 1/1028, when the 
quantization parameter r is changed from 0.0 to 3.0. The 
X 's represent the regular case, and show that the quanti- 
zation ambiguity there is large and energy-independent. 



The data points marked by squares correspond to the 
chaotic system, and clearly follow the linear law predicted 
by Eq. g (for d = 2). Finally, the plusses come from the 
"slow ergodic" sawtooth potential map ||l^ defined by 



V{q) = 0.3 



""2 



(13) 



where the quantization ambiguity is now put in the ki- 
netic term: H = ■ ■ ■ + rh^ cos 2Trp. Here we do not ex- 
pect quantum wavefunctions to be ergodic on the scale 
of a single channel in momentum space p5| , but they 
will be ergodic on scales ~ 1/| log?i| as ?i ^ 0. We then 
expect e ~ l/{a + P log h) in Eq. |^, which behavior is in- 
deed observed in Fig. 2. We see that ergodicity without 
chaos is sufficient to obtain unambiguous quantization of 
d = 2 systems, though the error in the high-energy limit 
approaches zero more slowly in the slow ergodic than in 
the fully hyperbolic case. 
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FIG. 2. Mean squared change in eigenvalues between two 
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as a function of the effective Planck's constant. 
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